Selective Mott transition and heavy fermions 
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Starting with an extended version of the Anderson lattice where the f-electrons are allowed a weak 
dispersion, we examine the possibility of a Mott localization of the f-electrons, for a finite value of 
the hybridization V. We study the fluctuations at the quantum critical point (QCP) where the 
f-electrons localize. We find they are in the same universality class as for the Kondo breakdown 
QCP, with the following notable features. The quantum critical regime sees the appearance of an 
additional energy scale separating two universality classes. In the low energy regime, the fluctuations 
are dominated by massless gauge modes, while in the intermediate energy regime, the fluctuations 
of the modulus of the order parameter are the most relevant ones. In the latter regime, electric 
transport simplifies drastically, leading to a quasi-linear resistivity in 3D and anomalous exponents 
lower than T in 2 D. This rather unique feature of the quantum critical regime enables us to make 
experimentally testable predictions. 
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I. INTRODUCTION 



Several years of intense experimental studies of quan- 
tum criticality in inter-metallic and heavy fermion com- 
pounds has lead to a growing evidence for a violation of 
the standard Landau Fermi liquid theory of the metal- 
lic behavior. Experimental phase diagrams use an ex- 
ternal tuning parameter, like the chemical pressure, the 
hydrostatic pressure or magnetic field to explore phase 
transitions at a temperature very close to the abso- 
lute zero. These phase transitions exhibit a regime of 
very strong quantum fluctuations -the quantum critical 
regime-, where anomalous thermodynamic and transport 
exponents are observed. A heavy fermion compound con- 
sists in a lattice of "big" atoms, like Ce , U or Yb, alloyed 
with a metal. The magnetic multiplet in the heavy atoms 
is determined by Hunds rules, and spin orbit interaction. 
For example, after spin-orbit interaction, the magnetic 
moment of Ce -4f 4 (S = 1/2, L = 3 ) is J = 5/2, for 
Yb-4 f 13 (S = 1/2, L = 3) we get J = 7/2 and for U 
5f (S = 1, L = 3 + 2) it leads to J = 4. The degener- 
acy is then split by the lattice crystal field effects, finally 
leading to a Kramers doublet. The apparent chemical 
complexity of these material has to be kept in mind for 
any investigation of their anomalous experimental prop- 
erties. In particular, the Anderson lattice model, which is 
standard in describing those compounds, relies upon the 
assumption that the Kramers doublets are well formed 1 - 
and thus, at low energy, one effectively considers that the 
physics is thoroughly described by f-electrons, sitting on 
the impurity atoms, subjected to strong Coulomb inter- 
actions , and hybridizing with the conduction electrons 
of the metal. The beauty of the experimental results 
lies in the fact that, although the compounds are based 
on a complex chemistry, very clean anomalous universal 
exponents in both transport and thermodynamics are ob- 
served. In this paper, we don't detail the experimental 



results, but rather refer the reader to the various reviews 
on the subject£i2^i. One of the most famous example 
of universality is the linear in T resistivity observed in 
YbRli2Si2 fine tuned to quantum criticality with a small 
magnetic field parallel to the c-axis. In this experiment , 
the linearity of the resistivity persists on three decades of 
energy, which makes it one of the most robust anomalous 
exponents in strongly correlated materials^. The coeffi- 
cient of the specific heat 7 = C/T is found to increase 
when getting closer to the QCP without showing any sign 
of saturation. A logarithmic law is generally attributed 
to this increase, but in several case, like for YbRli2Si2 
a divergence stronger than logarithmic is inferred from 
the data 3 -. Recently another type of experiment showing 
quantum criticality has appeared with the observation of 
strong effective mass increase in He 3 bi-layersi. In this 
experiments, one layer of highly compressed solid He 4 is 
adsorbed on a graphite substrate, then another layer of 
solid He 4 , and on top of it, two layers of He 3 are ad- 
sorbed. The first layer of He 3 undergoes a transition 
towards a localized state at a finite filling factor. Mea- 
surements of the specific heat show a strong increase of 
the effective mass inversely proportional to the coverage, 
associated with a decrease of the coherence temperature 
in (n — Jic) 1 ' 8 , as the QCP is reached. 

It is fair to say that this increasing body of fascinating 
experimental results is, at the present time, still myste- 
rious. The biggest challenge for the theoretician is to 
account for the linear resistivity in three decades of en- 
ergy, in 3 dimensions of space. Such a result calls for 
a specific scattering process, from which the lifetime of 
the conduction electron would acquire the characteristic 
linear in T dependence. At the moment, no theory can 
produce a reasonable scenario accounting for linear-in-T 
resistivity down to zero temperature, in three dimensions 
of space. The present paper, however, produces a very 
simple explanation for a resistivity linear in T , but in 
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FIG. 1: Schematic phase diagram for the selective Mott of 
the f-electrons in the Anderson lattice. On the right, when 
the holons condense, is the heavy Fermi liquid phase, also 
called in this paper the Higgs phase. On the left, where the 
holons are not condensed, is the localized phase for the f- 
electrons. We don't analyze here the magnetic character of 
this phase but rather characterize it only by the fact that the 
f-electrons localize. The QCP is multi-scale as is shown by 
the scale E* . Two regimes are of interest: (i) for T < E* the 
exponents depend on the shape of the Fermi surface, while 
(ii), for T > E* the exponents are universal and show quasi- 
linear resistivity in D = 3. To is the upper energy cut-off, 
above which the entropy Rln2 per site is released. 



the intermediate quantum critical regime; namely above 
a finite (although quite small) energy scale. 

Since the beginning of the history of heavy fermions, 
physicists have been struck by the intense magnetic na- 
ture of these materials. The presence of big atoms 
and big moments naively calls for magnetism. More- 
over, magnetic interactions mediated by conduction elec- 
trons - called RKKY interactions - are present in these 
materials^. They compete with the formation of the 
Kondo singlets responsible for the condensation of the 
heavy Fermi liquid phas o 12 ' 13 . For years the magnetic 
nature of the materials was so overwhelming that it was 
a real surprise when the first heavy fermion supercon- 
ductor was discovered!!. Similarly, the competition be- 
tween magnetism and the formation of the Kondo sin- 
glets has been considered until recently to be the main 
physical forces at play in the Anderson lattice 8 - 10 . In 
a completely natural way, the theory has first focused 
on QCPs separating a metallic magnetic phase from a 
metallic paramag net!!!£!£. This is this the simplest 
case of a QCP for itinerant electrons, usually called spin 
density wave (SDW) scenario. Here the Fermi surface 
of the conduction electrons is destabilized by scattering 
through magnetic modes propagating through the metal: 
the paramagnons. A key concept in this theory is the 
dynamical exponent z. Since close to a QCP quantum 
fluctuations are always relevant, the correlation in imag- 
inary time has to be taken into account which leads to 
an effective dimension d + z larger than the cZ-dimension 
of space. In practice, z is given by the structure of the 



bare paramagnon propagator. In the case of a transi- 
tion towards an anti-ferromagnet (AF), the paramagnon 
propagator has the form D~ 1 (q, iu n ) — c\uj n \ + q 2 where 
c is a constant. Hence z = 2 here. The theory then 
integrates the fermions out of the partition functio n 14 ' 15 
-step which is not rigorously justified -, to obtain an effec- 
tive <j) A lagrangian in d+ z dimensions. Since z = 2 in the 
AF case, one is typically above the upper critical dimen- 
sion of the model and the bosonic effective theory can be 
solved at the mean field level. A better treatment of the 
model doesn't integrate the fermions out of the partition 
function, but relics on an Eliashbcrg-typc theory where 
the vertices are neglected and self-energies conserve d 17 ' 19 . 
This theory is controlled in a rather artificial large N ap- 
proximation where N multiplies the number of fermions. 
However, we use it in this paper every time we deal with 
a QCP of itinerant electrons, since it is the best technique 
available at the moment in that it avoids integrating the 
fermions out of the partition function, hence not missing 
any infra-red (IR) divergence. Intense theoretical studies 
of the SDW scenario lead to the conclusion that some 
heavy fermion might fall into the universality class of the 
3D AF, like CeN^Ge^. Paramagnons, however, cannot 
account for strongly anomalous exponents like a linear 
resistivity in 3 D. 

The search for new QCPs then started. New ideas 
emerged, like the one of "de-confinement" of the heavy 
quasi-particle, close to the magnetic phase transition 3 . 
Meanwhile, a phenomenological approach was proposed, 
based on the observation that one can fit the NMR data 
for many compounds with two kinds of elementary exci- 
tations, that the authors called "two fluids"—. Although 
it's not clear whether it's really two "fluids" that are the 
good elementary excitations, the phenomenology looks 
robust and puts a landmark in the landscape of the field: 
the right microscopic theory should eventually reproduce 
the two fluids phenomenology. The first idea that a local 
mode appears at the QCP was given ir>22 in the so-called 
"locally quantum critical scenario" . The strength of this 
idea is to have been the first to outline that "something 
local" has to appear at the QCP. In this theory, the local 
mode is of magnetic nature (a local spin). Its weakness, 
however, is that it requires the magnetic structure of the 
material to be purely two-dimensional 23 . Moreover it 
appears to be very difficult to put it to the experimen- 
tal test, because the technique involved is not obviously 
transparent. To date, no experimental prediction has 
been made which would allow the community to accept, 
or discard it. 

Then, came the idea of a reconstruction of the Fermi 
surface at the quantum critical point, going with the idea 
of "de-confinement"—. The first tractable effective field 
theory illustrating the ideas of de-confincment and of 
Fermi surface reconstruction was given in 2 !. The authors 
used the Kondo- Heisenberg (K-H) model, and showed 
that, at some point in the phase diagram, there is a tran- 
sition towards a spin liquid phase, that they called FL*. 
At the transition, the Fermi surface re-constructs and 
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one looses half of the charge carriers in the phase FL*. 
Another way to look at the transition, is to say that, in 
the heavy Fermi liquid phase, the f-electrons start to con- 
duct, as was beautifully shown by 2 ^. Apparently the the- 
ory seemed to stand on much more solid grounds, since 
for the first time there was a tractable model on which 
to work. However, the treatment^ suffered some insuf- 
ficiencies. The most dramatic one is the complete lack 
of plausible experimental predictions. Apparently, the 
authors were not able to reproduce the linear in T resis- 
tivity, and the exponents they found for the anomalies in 
the specific heat coefficient didn't match the experiments. 
Moreover there was the problem of the phase FL*. To 
have a "spin liquid" phase in any heavy fermion com- 
pound appeared to be very unlikely. The most promis- 
ing candidate for it, was the system of He 3 bi-layers, for 
which intensive numerical studies in the last decad o 26 ! 27 
have shown with very little doubt that, in its localized 
phase, adsorbed He 3 layers form a spin liquid. 



II. RELATION TO PREVIOUS WORK AND 
STRUCTURE OF THE PAPER 

This paper details the theory of the "Kondo break- 
down", previously exposed in two letter o 28 ! 30 . We de- 
scribe the excitations around a selective Mott transition 
in the Anderson lattice, using the simplest theoretical 
tool available, the U(l) slave-boson gauge theory. Our 
goal is to cast the theory into the Eliashberg formalism 
with a careful treatment of the gauge invariance. To our 
knowledge, it is the first time that the Eliashberg tech- 
nique is presented for a "de-confined" QCP, where the 
order parameter breaks a gauge symmetry (instead of a 
global symmetry ). In particular, we present a study of 
the transport around such a QCP, using the Ioffe-Larkin 
(IL) composition rule. We derive as well the most generic 
Ward Identities (WI) associated with the gauge symme- 
try. 

First note that the first study of a phase transition 
separating a spin liquid from a heavy Fermi liquid phase 
was performed for the Kondo lattice model with frus- 
trated magnetism, using the DMFT technique^ 7 -. 

Two other studies exist. The authors of 2 ^ have shown 
that, at such a QCP in the K-H model, there is a jump 
in the transport properties; the residual resistivity jumps 
at the QCP and the Hall conductivity as well. This fol- 
lows from the fact that the f-electrons start abruptly to 
conduct in the heavy Fermi liquid phase. These results 
come naturally in our formalism, we thus agree with 25 . 
Our paper re-derives their results using the simpler for- 
malism of the IL composition rule. 

The second study is^ 4 -. The authors were the first to 
show that the volume of the Fermi surface jumps at the 
QCP. We agree with this result and with their observa- 
tion that the thermodynamics, in the vicinity of the QCP, 
is dominated by the gauge fields. We disagree with their 
findings about transport. The authors of ReL 2 ^ claim 



that at very low energy transport is dominated by the 
lifetime of the holons. One message of our study, is that 
transport is dominated by the conduction electrons in all 
regimes 2 ^. This is imposed by the IL composition rule 

P = Pc+ (pj 1 + Pb 1 ) 1 . (1) 

which states that the resistance of the holons is in series 
with the resistance of the spinons and in parallel with 
the resistance of the conduction electrons. 

Moreover, as already mentioned i n 28 ' 30 , the general 
structure of the QCP is much richer than what was found 
ir*21. This is depicted in Fi^TJ The first overlooked point 
is that the QCP is multi-scale. By this, we mean that an 
intermediate energy scale is present at the QCP, separat- 
ing two quantum critical regimes of universal exponents. 
The second overlooked point is that, in the intermediate 
energy regime, electrical conductivity simplifies, and one 
obtains a resistivity quasi-linear in T, in 3 D. This regime 
is of main importance, because it enables us to connect to 
experimental results and make predictions to test the the- 
ory. A first application of this regime to He 3 bi-layers has 
already been given 3 -^. The third overlooked point is that, 
at the mean field level, a modulated solution of the heavy 
Fermi liquid state exists, analogous in structure to the 
modulations found for FFLO superconductivity 32 . An 
analogous state was found in the study of Chromium 33 . 
This state may be of relevance for the mysterious phase 
observed in CeCoIns where the magnetic field is applied 
in the (ab) axis. Last, it appears that, in the case of 
the more physical model of the Anderson lattice with 
a small dispersion of the f-electrons, the Kondo bread- 
kdown QCP coincides with a Mott transition of the f- 
electrons^. 

The main advantage of this viewpoint is to simplify the 
discussion and to connect with other techniques avail- 
able, like Dynamical Mean Field Theory. DMFT studies 
can now be performed which will confirm or infirm in 
the short term the mean field findings of a selective Mott 
transition in the periodic Anderson lattice model (PAM). 
Selective Mott transition has been previously studied 
in the context of multi-orbital Hubbard model o 36 ' 37 i 38 
where transitions for various bands have been found. 
This model differs from the PAM because of the the ab- 
sence of hybridization between the bands. Recent DMFT 
studies 3 ^^^^. of the PAM tend to say that the selec- 
tive Mott transition indeed exist, but the studies are not 
completely conclusive yet. 

The second advantage is to link the discussion with the 
physics of high temperature superconductor o 42 ' 43 ! 44 ' 45 . 
It is notorious that the conduction electrons undergo a 
Mott transition in the phase diagram of high T c . Take 
a model with no frustration, like the K-H model on the 
square lattice. One can show formally that, in the lo- 
calized phase, the model is equivalent to a Heisenberg 
anti-ferromagnet and that its ground state is an anti- 
ferromagnet. This point was actually made in Re^ 4 -. Al- 
though this statement is formally simple, it is very dif- 
ficult to obtain within the slave boson technique. In 2 
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D, for high T c superconductors, one can advocate that 
the spin liquid naturally re-confines under the effect of 
gauge fluctuations^. We don't have this latitude in 3 D 
where gauge fluctuations are benign and can be treated 
within a non compact formalism. A better route towards 
the AF ground state is probably to allow the system to 
have both AF order and spin liquid, and then study the 
issue of re-confinement in that case. Similar studies have 
been performed in the early days of high T<4£ but have 
not yet been done for the Anderson lattice. 

On the theory side, the Kondo breakdown QCP suf- 
fers as well from its own weaknesses. We outline the 
main one, in our view, which is that this fixed point re- 
lies on the presence of a spin liquid at the transition. This 
spin liquid is described in terms of massless spinons, and 
the model is very sensitive to the Fermi surface of the 
spinons. To be precise, the exponents obtained using, for 
example, the uniform spin liquid arc different from the 
ones obtained using a nodal spin liquid. Whether the 
linear in T resistivity found in the intermediate phase 
survives the "spinonology" is still an open question. A 
related question is to ask whether the U(l) gauge the- 
ory used in the description of the Kondo breakdown is 
the adequate tool to describe the approach to a Mott 
transition. Are the elementary excitations correctly cap- 
tured within this technique ? This question is of broad 
interest and is as well at the heart of the physics of high 
temperature superconductors. 

The paper is organized as follows. In section IIHI we 
recall the model and introduce the slave-boson effective 
lagrangian. Section IIVI and [V] are devoted to the mean- 
field approximation. Precisely, in section [IV] we describe 
the mean-field approximation, showing the presence of 
the QCP at T = 0. In section [V] we recall the possi- 
bility of another mean-field solution: a modulated order 
parameter in the heavy Fermi liquid phase. In sections 
IVII to IVTiTl we study the fluctuations. We first start, in 
section IVII and IVII! by introducing the amplitude and 
gauge fluctuations within a "naive" RPA approximation. 
In section lVIIII we expose the Eliashberg formalism which 
links self-consistenly the amplitude and gauge fluctua- 
tions, and is the best available tool to study QCPs of 
that type. Precisely, in section [VI] we study the fluctu- 
ations of the amplitude of the order parameter, showing 
the different regimes of quantum criticality. In section 
IVIII wc introduce the gauge fluctuations, give the form of 
their propagator, and detail the gauging out of the the- 
ory. In section [Villi we recall the principles of the Eliash- 
berg theory and re-cast our study of the fluctuations in 
this framework. We then turn to electrical transport. In 
section IIXI we compute the Ioffe-Larkin (IL) composi- 
tion rule for the conductivity around the QCP and de- 
rive the various transport lifetimes in the low tempera- 
ture regime. In section [X] we focus on the intermediate 
regime and give the arguments for the quasi-linear resis- 
tivity in D = 3. In section 1X11 we give a summary of the 
thermodynamics and transport in this model. Finally, in 
section IXIII we present the conclusions and a criticism of 



the work. 

The Appendices are devoted to some technical details. 
Precisely, in Appendix [A] we give the calculation of the 
integrals used at the mean-field approximation. In ap- 
pendix (B[ we give the details of the evaluation of the 
vertices, which justify the Eliashberg treatment. In ap- 
pendix [D] we give an alternative, very simple derivation 
of the IL composition rules. In appendix [El we give the 
evaluation of the bosonic polarization, used in the form 
of the gauge propagator as well as in the holon trans- 
port lifetime. In appendix [F[ we give a field theoretic 
treatment of the Ward Identities of this gauge theory, re- 
lating any p-legs vertex to a (p-l)-legs one. We show as 
well in the most general way how the gauge symmetry 
protects the masses of the gauge fields in the Coulomb 
phase and how the mass is generated in the Higgs phase. 
In appendix [G] we give a direct check of the cancellation 
of the mass in the Coulomb phase, at the first order in 
the perturbation theory. Last, in Appendix [H] we give 
the derivation of ImS c and particularly show that the 
logarithm in D = 3 has a thermal origin. 

III. THE MODEL 

In order to study the Mott localization of the f- 
electrons in the Anderson lattice, we must first allow 
them to disperse. A small dispersion of the f-electrons is 
naturally present in the most physical models for heavy 
fermions. The reason why it is rarely included in the 
starting hamiltonians^ is that, in the study of the heavy 
Fermi liquid phase, the f-electrons dispersion is irrele- 
vant compared to the formation of the Kondo singlets, 
and one can safely approximate the narrow band by a 
flat one. We thus start with the Anderson lattice model, 
with a small dispersion of the f-band 

H = E ( C L^ JCT + fliaUj + /^,, •./•-) (2) 

where a is a small parameter, a is a spin index belonging 
to the SU(2) representation, Uj = t is the hopping term 
taken as a constant, V is the hybridization between the 
f- and c- bands, Eq is the energy level of the f-electrons. 

™f,i = YsafLfiv and n c,i = Y,a c L c 'i<? are tne opera- 
tors describing the particle number. We first study @ 
in the limit of very large on-site Coulomb repulsion U. 
U/ c accounts for the Coulomb interaction between the f 
and c electrons; although Uf c <C U it must be taken into 
account in the derivation of the effective theory. Here we 
go beyond the treatment o&. and include as well in the 
model the RKKY interactions mediated by the conduc- 
tion electrons. The starting hamiltonian writes 

H = E (4r«tf<*r + fUaUj + E S id )f ja ) (3) 



5 



E (V£U- + 



where Joi <A ar e determined by second order perturba- 
tion theory in large U/(at) and Uf c /(at) and Jrkky is 
determined by second order perturbation theory in small 
V/D where D is the bandwidth of the conduction elec- 
trons. One obtains J = 2(at) 2 /U, J\ — 2(at) 2 /Uf c 
and Jrkky = PoV 2 with po the density of states of 
the conduction electrons. S/ = J2 a /3 fa&apfp (resp. 
S c = c a&a0 c p) an d & is the Pauli matrix. The 

summation over (ij) can run on first nearest neighbors 
as well as second, third etc... nearest neighbors. Here we 
simplify the discussion by restricting all spin interactions 
to first nearest neighbors only, hence restricting our at- 
tention to the commensurate AF case. It has to be kept 
in mind that this is not the most generic situation and 
that frustration naturally appear in this model. 

We then take the U — ► oo limit of We account 
for the constraint of no double occupancy through a 
Coleman^ set of bosons (ir , b) enslaved to a constraint 

on each site J2a ftrficr + b\bi = 1- We make at each site 
the transformation 



fi<7 * fi 



rb] 



(4) 



where the f*- creation operators are called "spinons" and 
the "holons". The physical electron / splits into a 
"holon" and a "spinon" under the transformation ([4]). 
We notice that, upon this transformation, the slave boson 
drops out of all bi-linear products of fields at the same 
site. Indeed, one can explicitly show that at site "i" : 



flhb\U\0) = flfia\0) = 
flWfiM = flUW) = n f , ia \a) 



(5) 



Another way to see this is to apply the constraint inside 
the operator. 



= (1 - n b )flj ia 



(6) 



The last line is obtained by observing that the constraint 
of no double occupation is equivalent to n 2 ^ = nfi. 
The effective Lagrangian is then 

£ = E {4a( d r Sij +Uj) Cja (7) 

+ CAI>,"I>jI>j + (dr + E + Ai)^)/^) 



{ij) 



where the constraint has been implemented through a 
Lagrange multiplier A^. Note that the spin operator 
S/ = 2 a/3 fa&apfp is now expressed solely in terms of 
the spinons and thus is insensitive to the slave bosons. 

In the following of the paper we consider a "large N" 
extension of the Lagrangian |(7J) by enlarging the spin 
group from SU(2) to SU(N). The indices a now belong 
to the SU(N) group. 



IV. MEAN-FIELD APPROXIMATION 

In the mean-field approximation, we minimize the ef- 
fective action with four fields : (A, b, <fro, erg ), where a 
static and uniform approximation is made on all fields. 
4>q is the uniform spin liquid parameter, which de- 
couples the short range AF interaction JS ■ S — > 

( t ) oYji,a(fiafi.<y + h - c -) ~ N( t>o/ J - a Q is tne uniform 
field which decouples the induced Kondo-like interaction 

JiS/, v S,., - aoEUfl^ + h.c.) - Nal/ A. The 
minimization of the free energy leads to the following 
mean- field equations: 

T E ae k G„ + j E G fc + e f - E = ; 

T E G fc + Na /Ji = 0; 

T E (e k /D)G ff +Ncf> /J = ; 

b 2 +T E G // = ^7 2 > 



where e k is a typical dispersion of the conduction elec- 
trons, e/ = Eq + A and the dispersion of the spinon band 
is taken to be e° = (ab 2 + (4>o/D)) e k + £/■ We have the 
following Green's functions 



G 



ff 



'fc 



G c 





lLU n 


- £fc 




(iLU n 


- e k )(iuj n - 


4) - 


(Vb + ao) 2 




Vb- 


+- c 




(iuj n 




4)- 


(Vb + aa) 2 






-4 




(iu> n 


- e k ){iiv n - 


4)- 


(Vb + a Q ) 2 



(8) 



The mean field equations are solved in the case of a lin- 
earized bandwidth and for N = 2; energy scales extracted 
from the mean -field studies are as well written for N=2. 
The summations over (k,uj n ) can then be performed an- 
alytically and are given in Appendix \Kl the set of re- 
sulting equations is then solved numerically. Since the 
interaction JiSf^ ■ S CJ - generates some additional Kondo 
coupling, it is not obvious that a Mott transition- a point 
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in the phase diagram where 6 = 0- occurs at finite V. 
However this is what happens, the additional Kondo cou- 
pling Co being itself driven to zero at the transition. The 
result 
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FIG. 2: Effective hybridization Vb (black curve), the f-band 
chemical potential e/ = Eq + A (blue curve) and 2O0o (red 
curve), as a function of V. The electron bandwidth is D — 
1000. The chemical potential fi = 0, the ratio of f- and c- 
masses is a = 0.1. <j>o is evaluated self-consistently around 
finite RKKY value of 10" 3 D. The f-energy level is E = 
— 500. The mean field equations are solved for N = 2. The 
effect of the field a is negligible, hence not shown in the figure. 

We can study analytically the nature of the the fixed 
point by expanding the mean-field equations {8]| around 



cto = ; 
6 = ; 

1 0o I = cJ; with c ~ 0.1 . 
From the first equation, we get that 
n /c (o) + (e f - E )/V 2 ~ 

1 



(9) 



where U fc (0) = T ^ 



(10) 



k,0J n ,<7 



(iw n - e k )(iuj n - e9) 



Using, for the conduction electrons and the spinons, 
linearized dispersions of the form e k = vp(k — k F ) , 
e 1 = v o(k — ko) with vq = 24>o/kp one obtains 



with a' = ab 2 + <fi /D 



(11) 



(12) 



One can convince oneself that the fixed point occurs when 



= DExp 



Eg 

N P0 V 2 _ 

Kondo scale of the problem 50,5 ! . Hence 



One recognizes here the typical 



Jcrit = -DExp 



En 



Np V 2 



(13) 



One sees that the location of the QCP in the phase dia- 
gram depends crucially on the spin liquid parameter J. 



Moreover, it is important for the stability of the mean 
-field equations that the bandwidth stays finite when 
6^0. This is achieved with the uniform spin liquid 
parameter (j)Q being finite at the transition. The mean 
field equations are of little help to determine the nature 
of the localized phase where 6 = 0. In the simplified 
scheme we have taken, this phase is a uniform spin liquid 
- 4>o being finite and uniform. The spinons can however 
have some other symmetries, and it would be useful, for 
example, to determine the location of the transition when 
spinons with a nodal Fermi surface are used. One can as 
well allow for AF order in the localized side and study 
the stability of the phase diagram. This is the program 
for future work. 



MODULATIONS OF THE ORDER 
PARAMETER 



In the previous section,we have naturally considered 
that the order parameter b was ordering at q — 0. This 
issue has to be re-considered keeping in mind the nature 
of the spinon Fermi surface. In this section we review for 
completeness, the results obtained in2£. At the QCP, the 
effective mass of the order parameter writes 



D-\q,0) = m b ; 

m b (q) = Po (-E + V 2 Il fc (q)) 



(14) 



where Hf c (q) is the static /c-polarization, taken at fi- 
nite momentum q but zero frequency. At the QCP , the 
minimum of the effective mass determines the ordering 
wave vector. Two situations are to be considered. First, 
if both the f- and c- bands are electron-like, Hf c (q) ob- 
tains its minimum at qo = 0. Note that in that case, the 
curvature of the Fermi surface has to be included to see 
the minimum. Second, if one band is electron-like but 
the other one is hole-like, then the situation is analogous 
to the FFLO ordering in superconductors^ or to change 
density waves in Chromium 3 -. At the QCP, the ordering 
wave vector is at go — 1.2(7* where 



q* = \k F - fc | 



(15) 



is the difference between the Fermi wave vectors of the 
two species. The determination of the ordering wave 
vector inside the ordered phase is much more involved 
than at the QCP. It led to a full literature in the case 
of FFLO superconductivity^. One generally expects a 
first order transition towards a uniform order inside the 
ordered phase. Here the situation is rather more com- 
plex than in the FFLO case because the order parameter 
6 carries one quantum of gauge charge. Gauging out the 
theory, even at the mean field level, is required to deduce 
the observable quantities. This study definitely deserves 
more work. Particularly, it would be interesting to see 
what kind of superconductivity occurs in a Kondo phase 
where the hybridization is modulated in space. 
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VI. AMPLITUDE FLUCTUATIONS 

In the two following sections, we describe the random 
phase approximation (RPA) evaluation of the amplitude 
and gauge fluctuations. The fluctuations are studied in 
the case of the order parameter condensing at q = 0. The 
fluctuations of the amplitude of the order parameter are 
more complex than what was considered ir>21 since the 
order parameter is coupled to two types of fermions: the 
f-spinons and the conduction electrons. Neglecting for 
the moment the effect of gauge fluctuations, within the 
RPA the polarization is similar to a Lindhard function, 
but with two different types of fermions. 



-Po 



(17) 



2vp a'q(l — a') 
[(-a'iCl n + a'v F {q - q*))Log(-ia'Vl n + a'v F (q - q*)) 

— (— a'if2 n + a'v F (— q — q*))Log(—ia'Cl n + a'v F (—q — q*)) 

- (-iCl n + a'v F (q - q*))Log(-iCl n + a'v F (q - q*)) 

+ (—ifl n + a'v F (-q - q*))Log(-itt n + a'v F {-q - q*))] . 

We see that an additional scale q* naturally emerges from 
the polarization, q* is the difference between the Fermi 
wave vectors of the f-spinons and c-electrons. Here we 
have taken both the spinon Fermi surface and the elec- 
tron Fermi surface to be spherical and centered with re- 
spect to each other, hence q* isotropic. We define the 
energy scale for N = 2 by 



n 



fc 



E* = 0.1a' D 



k F 



(18) 



We can expand (fT7|) in four different regimes, (i) For 
q < q*, |O n | < E* we have 



D b l (q, iCl n ) = po [—iCl n + S + aq 2 + H/ C (q, ifi n )] , 

(16) 

with a — — iVLog(a')/[(l — a') 2 k F ], the bosonic mass 
<5 = — Eq comes from the mean-field equations, and p is 
the density of states of the conduction electrons. Note 
that, in our definition, If/ C is the "dynamical" polariza- 
tion, corresponding to the infra-red (IR) sector. In the 
text, we use this definition for every polarization. Let's 
here evaluate II/ C . There are two cases of interest, as 
shown in Figj3] Depending whether the Fermi surface of 
the spinons and conduction electrons intersect or not, we 
have two forms for the amplitude propagator. Let's start 





Case a) 



Case b) 



FIG. 3: Illustration of the two typical cases of interest. In 
case a) the two Fermi surfaces of the spinons and conduction 
electrons are centered and there is a gap between them; in 
case b) the two Fermi surfaces intersect. 

with case a) where there is a gap between the spinon 
and electron Fermi surfaces. Using linearized bands, the 
polarization can be computed analytically and we get at 
T = 0: 



Uf c (q,i£l n ) 



Npo 

47T 



de duj dcosO x 



(iuj + ifi n — e — v F q cos9)(iuj — ae — av F q*) ' 



D b 1 (q,iQ n ) ~ po 



~ 5 + a( f- N J^_ 
a'v F q* 



(19) 



with 5 = -Eo + p V 2 Logol / '(I 
q* , E* < |fi„| we have 



a'), (ii) For q < 



D b 1 (q,ifl r , 



Po 



aq2 + N Lo g \n n \ 

a'v F q* 



(iii) For q>q*, a'v F q < \Q n \ < v F q we have 



D b 1 (q,iil n ) ~ po 



aq2 + N Log\n n \ 

a v F q 



(iv ) For q > q* , v F q < |0„| we have 



D b 1 (q,ifl rl ) ~ po 



aq 



N- 



10, 



a'v F q 



(20) 



(21) 



(22) 



Note the peculiar form of the boson propagator, where 
the different sectors are not in the same footing in large 
N§&. We don't know any way to avoid this problem 
within the Eliashberg theory, and that's why we call it 
a large N expansion rather than a large N limit. We re- 
fer the reader to 19 ' 56 for further details about this and 
to next section, just mentioning here that the merit of 
the Eliashberg theory is to capture all the frequency di- 
vergences within the leading and second leading orders in 
1 /N. To each regime one can associate a different dynam- 
ical exponent. In regime (i) the dynamical exponent is 
z = 2 and the propagator corresponds to an un-damped 
bosonic mode. In regime (ii) the dynamical is exponent 
z = oo. In regime (iii) we get z — 1. Finally, regime 
(iv) has z = 3. It was shown that the spectral weight in 
the (f2, g)-space is most entirely centered in the regime 
(iv) where z — 3 28 . This feature is illustrated in Fig 
21 Physical arguments enable to simplify the discussion. 
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momentum q 



FIG. 4: Spectral weight of the n/ c polarization in the (q,uj)- 
plane. We have taken q* = 1, vf = 1, a 1 = 0.01. We see that 
the weight is concentrated in the regime II (z — 3) for q > q* , 
where one has particle-hole damping. 



At low momentum and low energy, the particle-hole con- 
tinum is gapped, due to the mismatch of the two Fermi 
surfaces. There is propagation of a single boson mode 
with exponent z — 2 (|19p . At high momentum and high 
energy, one can neglect the mistmatch of the two Fermi 
surfaces; the two fermion species behave as if they were 
identical. The polarization behaves like a Lindhard func- 
tion ([22]) with z = 3; this form is typical of a q = phase 
transition. 

As we have shown above, an energy scale E* is present 
in the quantum critial regime, separating two regimes 
with universal exponents. In the low energy regimes 
T <C E* the fluctuations of the amplitude are gapped 
and thermodynamics is dominated by gauge fluctuations. 
For T 3> E* the thermodynamics is in the universality 
class of the ferromagnet with z = 3. 

On can ask the question: -"why is E* so low ?" To 
answer it we give an estimate of E* . One can first iden- 
tify a' D with the temperature To at which the entropy 
i?Ln2 is released in the Anderson lattice. This scale 
is seen experimentally in every compound, as a bump 
into the thermodynamic and transport observables. For 
T > a'D, the spinons loose their dispersion; hence one 
can consider that they behave as free spins and that the 
entropy RLn2 is released. This observation fixes the scale 
To to be roughly 24K for a compound like YbRIi2Si2 and 
50K for CeCoIn 5 . We see from Eqn. (fl"8)) that even for 
very large q* , E* is already an order of magnitude lower 
than To, which makes it a small scale already. Next, de- 
pending on the form of the spinon Fermi surface, q* can 
be as small as O.lfcF, in which case the scale E* is of 
the order of lOmK; which means it is not experimentally 
reachable. The fact that E* relies on the imponderable 
form of the spinon Fermi surface is assuredly the weak 
point of this theory. 

To illustrate this point, we look at the more generic 
case b), where the spinon Fermi surface and the elec- 
tron Fermi surface intersect in two "hot lines" . Since the 



anisotropy is broken, we define a new q* as 

q* = Max (M) p F -fco|] , (23) 

where the maximum is taken over anglular variables on 
the Fermi surface. In this case, the regime (iv) is un- 
changed, but the regime (i) is different, since the particle- 
hole continum is not gapped. The boson propagator in 
regime (i) is now 



D b 1 (q,in n ) ~ p Q 



aq 



Nc- 



a'vpqo 



(24) 



where qo is the modulus of the wave vector at which the 
two Fermi surfaces intersect, and c is a cocfficcnt depend- 
ing on the shape of the Fermi surfaces. We see that the 
propagator still has the dynamical exponent z = 2 but it 
is now damped; hence in the same universality class as 
a 3 D anti-ferromagnetic SDW QCP. We note that the 
change of the spinon Fermi surface between case a) and 
case b) didn't affect the regime for T > E* . Although 
we cannot prove that this result is general, we illustrate 
it with considering a third case where the f-spinon Fermi 
surface is nodal in 2 D, as shown in Figj5] This case is 
interesting because it is the exact analogue to the descrip- 
tion of the one band Hubbard model within the U(l) or 
SU(2) slave-boson gauge theor y 42 ' 45 . One can show that 
the amplitude propagator is formally equivalent to case 
a), where the two Fermi surfaces are gapped. Hence in 
that case as well, the regime with T > E* is unaffected 
by the form of the spinon Fermi surface, although E* is 
supposedly bigger that in case a), since q* and /cf are of 
the same order of magnitude. In conclusion, the regime 
T > E* depends on the existence of the spinon Fermi 
surface but not on its shape. 




Case c) 



FIG. 5: The third case of interest, where the spinon f Fermi 
surface is nodal, like for the U(l) gauge theory of the one band 
Hubbard model, in D — 2. This is in the same universality 
class as the case a) 



VII. GAUGE FLUCTUATIONS 

The order parameter of our phase transition is not a 
macroscopic order parameter, but relies on the conden- 
sation of a gauge boson carrying one quantum of charge 
of the gauge invariant theory. We study here the U(l) 
local invariance of the Lagrangian. For clarity we start 
with the evaluation of the polarization of the gauge fields 
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within RPA. The WI constraining the masses are given in 
Appendix[F] We start with the version of §5§ after having 
performed the Hubbard-Stratonovich transformation. 



"3<r 



C = ( c L( d T 5 lJ +t lj )c 3 

+ flikatijb] + {d T + E a + iXi)Sij)fj. 

i 

N, „ 



Ji 



Under the gauge transformation 

fi * fi e 1 ? 
bi^bi e ie ;, 
Oi -> Oi e l9 ' , 
A; — > Xi + d T 6i, 
6,, -> 6u e* 9 -- 9 - 



(25) 



(26) 



the Lagrangian (|25|) acquires a total derivative C(t) — > 
£(r) — i^iidrOi) which is a multiple of 2i7r and thus 
leaves the Lagrangian invariant. The above Lagrangian 
is thus invariant under a U(l) local -or gauge - symme- 
try. In the mean- field treatment above, we have consid- 
ered that the Hubbard-Stratonovich fields can be taken 
at their saddle points (cr*) = <to, (Xi) = iX , (6j) = b 
and (<pij) = 4>o and we have neglected the fluctuations of 
the phases of the field as well as of the amplitudes. It 
is safe to assume that the amplitude fluctuations of <f>ij 
are gapped since 4>q doesn't vanish through the phase 
diagram. At the QCP, however, the amplitude fluctu- 
ations of bi and ai becomes massless. Another issue is 
the phase of the gauge fields. In order to study them, 
it is convenient to take the continuum limit in the La- 
grangian (|25p . For this, we write <pij — 0oe laii and now 
the second and third terms in Eqn. (|25p describe spinons 
and holons in the presence of the fictitious electromag- 
netic field B associated with the vector potential a with 
J 3 a dl = SLi j 42 ' 45 . We then coarse grain the Lagrangian 
to obtain the continuous limit 



£(r,r) = E/ d ^l{dr 



2m 



A* c 



t , (V + ie a/c) 2 , 

"° ' 2^ f m ° I /<T 

(V + ie a/c) 2 



6 f (dr 

N > ,2 N 



2mb 
'I 



+ X + ia ) b 



£ J dv[{Vb + a)fic„ + h.c] 



(27) 



When the order parameter b + a/V condenses, the gauge 
field a acquire a gap proportional to (b + cr/V) accord- 
ing to the rules for the condensation of a Higgs bosor*^. 
This phenomenon is responsible for the Meissner effect 
in the case of a superconductor. In Appendix [F] we de- 
rive the generation of the mass using the Ward identities 
associated with the gauge invariance. 

A non trivial issue is the gauging out of the theory. The 
slave-boson representation required the use of additional 
fields, hence the field theory is redundant and some gauge 
fluctuations can be factorized out of the partition func- 
tion. Since one deals with a U(l) gauge theory, one gauge 
fixing constraint only is required. When gauging out the 
theory, it is convenient to chose the radial gauge for the 
Higgs phase - also called here "heavy Fermi liquid" , which 
follows previous studies of the Kondo impurit y 48 ' 49 and 
of the Kondo lattice 1 ^. Namely, we choose the order pa- 
rameter to be real, with for example 



b + a/V = \b + a/V\ 



(28) 



This constraint is alternatively called the "physical 
gauge" in field theory, since the gapping of the gauge 
fields in the Higgs phase is transparent in this gauged. 
When the order parameter vanishes, in the so called 
Coulomb phase, it is clever to use the Coulomb gauge 
imposing the condition 



0, 



(29) 



so that the vector fields become purely transverse. In the 
Coulomb gauge, the scalar fields (/i = 0) and vector fields 
decouple; the scalar fields are massive, which are nothing 
than the density-density correlation function, while the 
mass of the vector fields remain massless (see Appendix 
IF)) . The polarization of the vector fields gives rise to a 
Reizer singularity^ described below. Of course, formally 
all gauge fixings should be equivalent, and there is no 
deep reason to chose one constraint rather than the other. 
However, if for example one tries to work with the radial 
gauge in the Coulomb phase, lots of divergences appear 
which have to cancel at the end, for physical quantities. 
The reason is that radial coordinates are ill-defined when 
the radius vanishes (here the radius is the modulus of the 
order parameter |6 + cr/U| ) . Some divergences are then 
hidden in the Jacobian of the transformation. We are 
not aware of any field theoretic treatment of this case for 
condensed matter systems. 

At the QCP we thus work in the Coulomb gauge where 
the vector fields have the purely transverse form : 

D~Mq,iQ, n ) = Tl(q,iCl n )(8ij - qiq 3 /q 2 ) ,with (30) 



j(T T [J / , i (r,r)J /)J -(0 > Q)]) 



10 



2m f 



-(T T [J M (r,r)J bii (0,0)]) 



+ & 



Po 
2m b 



5(r)S(r) 



(31) 



The first term in (|31[) corresponds to the paramag- 
netic contribution while the second term is the diamag- 
netic part. J/ is the current operator defined as J/ = 
^/(2mo)[/ t V/~(V/t)/] and J b = i/(2m b )[tfVb-{W)b] 
After Fourier transforming, the polarization can be writ- 
ten as 



Ii i3 {q,iVL n ) 



^2 E (k+q/2)i(k + q/2) j 

P f 

x G//(k + q, iw n + ifi„)G//(k, iw n ) +%^^ 
T 



—^(k + q/^k + q^ 



x G b (k + q, iw„ + i£l„)Gb(k, iuj n ) + % 



2m b 



The vector field propagator has the form 

D^{q,m n ) = (11/ + n„) - <m,Az 2 ) 
where at the QCP and at T = 0, 



(32) 



II/(q, iO„) 
n & (q,ifi n ) 



iV 
2m/ 
1 

2m b 



{q/k F f 



(33) 



with / d = J( ( ; 2 /d)d( d - 1 ) ( 7/ (q)/(2^)^ 1 and / (g) is a UV 
cut-ofS^. The evaluation of lib is given in Appendix [El 
Pictorially, the polarization can be written 



/ 



II/ + life — w» <tw 




VIII. ELIASHBERG THEORY 

Our QCP corresponds to a q = transition. In the 
two previous sections we have treated the fluctuations 
within the RPA, for which the polarization is evaluated 
at the first loop order. A self consistent treatment of the 
fluctuations is needed to control the results. Integrating 
the fermions out of the partition function is a dangerous 
uncontrolled step for such transitions. A better approach 
is the Eliashberg theory, controlled in a large N expan- 
sion. For the details we refer the reader to the exten- 
sive review of this technique given ir>i2. For complete- 
ness we recall here the reasoning and the main results. 
The Eliashberg theory relies on three steps. The first 
step is to neglect the renormalization of the vertices as 
well as the momentum dependence of the self-energy. In 
the second step, the Dyson equation is used to evaluate 
self-consistently the boson polarization and the fermion's 
self-energy. Then one checks that the approximation is 
correct. Here we have two types of fermions and two type 
of massless bosons as well - the order parameter and the 
vector gauge fields; the time gauge field ao being massive. 
The coupled Dyson equations are represented below with 
Feynman diagrams 

Gj 1 (k,uj n ) = iu n - eo + iS/(w n ); (34) 
G~ 1 (k,uj n ) = iu n - e k + ffi c (w„); 

D b 1 (q,^l n ) = Po [-i&n + 5 + aq 2 + U fc (q, O n ) - S 6 (O n )] ; 

D^( q ,n n ) = Po [{q/k F ) 2 + n f { q ,n n ) + ii b {q,n n )] 

x (S i:j - q t q /q 2 ) , 



where 



S/(w„) 



£ c (w„) 



f c 



b, a M 



[TRnnr] 1 = [wvv] 1 + [vwv] ; 



Sb(^n) 



When the bosons condense the gauge fields propagators 
become massive. Note that the two-legs vertices in a M - 
are protecting the mass sector—. This is shown in 
Appendix [Fl using again the Ward identity associated to 
the gauge invariance of the problem. It is checked as well 
in Appendix [G] by a direct evaluation at the first order. 



II/ c (<j, Q n ) 



/ 
/ 

n/(g,n n ) = 
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b 



U b (q, n n ) = 




b 



Note that, in the above diagrams, the lines are full by 
construction, but the evaluation of the diagram doesn't 
change if the fermion propagators are taken to be bare. 
The bosonic "self -energy" £5 and the polarization lit, 
are given in Appendix [E] 

The crucial point in the Eliashberg theory, is to check 
that the three-legs vertices are small. We recall that, 
in the regime T > E* , the dynamical exponent is z — 
3. This exponent characterizes as well the ferromagnetic 
QCP 19 and the U(l) gauge theories 06 , which have been 
studied in the literature. In both cases, the Eliashberg 
theory is controlled in the same way, using the combined 
effects of the large N expansion and the curvature of the 
Fermi surface. We recall here the results. 

There are two types of vertices depending on the in- 
coming momentum and frequency. The static vertex, 



T(0,0) = 




k F ,0 

obtained for vanishing incoming frequency, is small in the 
large N limit, going like T(0, 0) ~ iV -1 / 2 . The dynamical 
vertex, however, obtained in the limit of non vanishing, 
but still small frequency when \Q n \ <C t'F<7 is dangerous. 
Note that within our model it is impossible to form the 
mixed /c-vertex at the first order. Hence we form it at 
the second order, but evaluate each loop separately in 
Appendix [B] 



r(g,fi) = 




k F ,0 

For bare fermion legs, it behaves as T(q, O n ) ~ Sli 3 d ^ 3 
(for d = 3 we get T(q, f2 n ) ~ — fi n Log|fi n | ) hence show- 
ing strong anomalies in the frequency dependence in the 
infra-red sector. Fortunately, in the Eliashberg theory, 
the fermion legs are not bare, but dressed with their 
proper self-energy. For z — 3 the spinon and electron self- 
energy behaves as £(w n ) ~ \u; n \( d ~ 3 ^ 3 (with for d = 3 , 
£(w„) ~ — w n Log|f2 n | ). Dressing the fermion legs effec- 
tively kills the divergence in frequency. A subtle point is 



that, with a linearized Fermi surface, the resulting ver- 
tex is of order one, with no small parameter. Includ- 
ing the curvature of the Fermi surface changes this sit- 
uation. In the case where qy/m* 3> Vfq x where m* is 
the curvature mass, we finally get T{q, 0„) ~ (3 2 LogP 2 
with (3m*/(Nm) <C 1. As we see, the justification of 
the Eliashberg theory in the case ofaz = 3,t7 = 
QCP is not simple; it requires both the presence of the 
curvature and an additional large N expansion (N being 
here the number of species of fcrmions), which is rather 
artificial. An additional caveat is that the large N limit 
cannot be rigorously taken since, if we did this, we would 
not include the diagrams corresponding to the fermionic 
self-energy which scale like TV -1 / 3 in D = 2, hence are 
sub-leading in the large N limit. The Eliashberg theory 
thus relies on a large N "expansion" rather than on a 
"large N limit" . This remark follows the similar observa- 
tion made in the previous section, that the form of the 
boson propagator is not homogeneous in N. 

Another regime of interest is the case where z = 2 , for 
intersecting Fermi surfaces. The theory is equivalent to 
the AF QCP whose Eliashberg spin fermion treatment is 
reviewed inii. We refer the reader to this paper to get 
convinced of the smallness of the vertices. 

The two regimes (ii) and (iii) of the boson propaga- 
tor are not physically relevant since the boson propaga- 
tor has only a small spectral weight in this regime com- 
pared to regime (iv). We still must check that the ver- 
tices are small. Let's check the regime (ii) where the 
frequency part has a characteristic logarithmic logarith- 
mic dependence. In both 2 D and 3D we find that the 
static vertex scales like F(0,0) ~ \/N. Curvature is 
needed to regularized the dynamical vertex, which goes 
like T(q,Q n ) ~ (f2 n Log|Log|f2 n ||)/Log|f2 n |. Hence large 
N is not needed for the smallness of the dynamical vertex 
in this regime. Details of the vertex calculation are given 
in Appendix [Bl 

Last, the reader may wonder whether the analogous 
of the singularities discovered by Belitz, Kirkpatrick and 
Vojta (BKV}i£ exist in this theory. This issue lies be- 
yond the Eliashberg theory, but is of importance for the 
stability of q = QCPs. It has been shown by the au- 
thors of 18 that singularities appear in the static sector, 
in the computation of the static, and temperature de- 
pendent polarization. Indeed close to a q — QCP of 
ferromagnetic type, the static polarization 



c 



I%T) = 




c 



vanishes at the first loop order. Technically, one can 
check that, with zero incoming frequencies, poles of the 
fermion lines are in the same half-plane; hence the static 
polarization vanishes at one loop. To get a non zero result 
one must include the first vertex correction. 
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In our case, it is worthwhile to notice first, that the 
one loop polarization has a temperature dependence at 
lower energies, due to the fact that we have two fermion 
species. Indeed we find by direct computation an acti- 
vated behavior in the case of gaped Fermi surfaces: 



n /c (T) = n 



(1-aO 



Up 



a vpq 



(1-a') 



where np is the Fermi distribution function. This acti- 
vated behavior is quite small. We have as well a source 
of damping coming from the gauge fields and leading to 
the temperature dependence (see AppendbiE| 

E 6 (T) ~ T< d+2 )/ 2 . 

However the BKV singularity comes from inserting ver- 
tices and self-energy (beyond Eliashberg theory) in the 
polarization bubble. It is dangerous because the first in- 
sertion leads to a contribution with negative sign, hence 
de-stabilizing the QCP. Here we first observe that there 
is no first order vertex correction to the polarization Hf c ', 
the diagram simply cannot be formed. At the first order 
one has only self energy insertions of the type 



ng(T) 



b /"WvA b 



b /^V\W\ b 



We find, in the intermediate regime, in D = 3 

n W ~ _ N -^/3 T 4/3 

We argue that this contribution is smaller than E* and 
can be neglected. In the low energy regime we find that 
the mismatch between the two Fermi surfaces protects 
the system, independently of whether the Fermi surfaces 
intersect or not. 

The first vertex correction comes at the two loop order 
as depicted below. Note that, contrarily to the nematic 
case^, there is no cancellation between the vertices and 
self energy insertions. This is due to the fact that we have 
two fermion species, hence the cancellation is spoiled, 
even though the transition is in the charge channel. We 
find that this correction behaves as 



in the intermediate regime in D = 3, and has a positive 
sign. It is thus not dangerous. 



nf c {T) 




(Na')- 2 / 3 T 5 / 3 



Note that the fermion lines are full in the computation 
of these diagrams. The interested reader can find more 
details about this discussion in Appendix [C] where we 
comment as well on D = 2. 



IX. IOFFE-LARKIN (IL) COMPOSITION 
RULES 



In this section we start the study of transport prop- 
erties of this quantum critical gauge theory. A previ- 
ous study exists 2 ^, but here we re-cast the formalism in 
terms of the Ioffc-Larkin composition rules for the resis- 
tivity and enlarge the discussion to the temperature de- 
pendence of the resistivity . It is known from the seminal 
paper of Ioffe and Larking that in the Coulomb phase 
of a gauge theory, both spinons and holons participate 
to transport and the total conductivity is constrained by 
the Ioffe-Larkin (IL) composition rules. In order to de- 
rive them for our model, we must expose the system to 
an external electro-magnetic field A. A is attached to 
the conduction electrons and to the holons b^. The sys- 
tem is as well subjected to the internal fictitious gauge 
field a (f2"T|) . a is attached to the fields carrying the gauge 
charge, hence to the holon b an the spinon We work 
in the Coulomb gauge, ensuring that V • (A + a) = 0. 
The system is invariant through: 



A -> A + V6 A 



a -> a + V6> a 



c -» c Exp(i0 A ) / -» / Exp(i6 a ) 

b-fb Exp(i£» A - i8 a ) . 



(35) 



In order to derive the IL composition rules, we expand the 
action for minimal coupling with the vector gauge fields. 
Note that the composition rules are completely general 
(see Appendix [Dj ; the expansion in minimal coupling is 
a matter of convenience. 

S = f drJ^PGj 1 ^)/ (36) 

J k 

- I dT E G »-J f + ct °o.c c + & f ^oj b 
J k 

+ p v/ • a / + c f v c • A c + 6 f v b • (a + A) 6] , 

where the vertices v/, v c , v b are determined from the 
expansion of the effective lagrangian to the first order in 

a, A£2. After integrating out the matter fields f,c and 

b, the effective action for minimal coupling to the gauge 
fields reads: 

S{A,a} = ^J ^^[A Q (c,fc)H^( w ,fc)A /3 



(2n) a 

(A Q -a a )nf (w,fe)(A /9 -a fl ) 
a. a (ui, fe)lK (a;, k)ap 



(37) 



The polarization bubble H(cj, k) can be decomposed into 
its longitudinal and transverse parts: 



n a /3 = S a 



k a k^ 



(38) 
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Following IL, we integrate over the fictitious gauge field 
a in © to get the effective action for the external electric 
field. 



S{A,a} = j I dkJ2^ a (u,k) Ile + ^+nj 1 )- 1 



Hence the total polarizability of the system is 



p l = n c + ( ny 1 + nr 1 



Using the Kubo formula 



(39) 



(40) 



(41) 



we get for the conductivity 

a = (Jc + {vj 1 + O" 1 • ( 42 ) 

Note that the composition rule for the conductivity can 
be derived in a simpler way, using solely the constraint 
attached to the gauge symmetry; this calculation is given 
in Appendix [D] The composition rules (|4"0")) are power- 
ful enough to allow a complete discussion of the electrical 
transport close to the QCP. One observes that the holons 
are "sandwiched" between the spinons and the conduc- 
tion electrons. If 11^ 1 dominates II J 1 , then it is very 
unlikely that Iff, will dominate II C since II C already dom- 
inates IT/. Hence we infer from the simple form of (j4"0|) 
that n c is the most important part. The transport is 
dominated the conduction electrons. 

Let's first examine the limit of zero temperature. On 
the heavy Fermi liquid side, the holons b are perfect con- 
ductors and the spinons are massless fermions so that 



H, 



n 



/ 



n/e 



(43) 



mt rrif (1 — iu/Tf) 

where we have used a Drude formula^ for the polariza- 
tion of the f-spinons and r/ is the scattering lifetime of 
the spinons^, which writes tJ 1 ~ Tq 1 + T 2 in the heavy 
Fermi liquid. Note that we have taken into account the 
effect of impurities in the scattering time of the spinons. 
This corresponds to dressing the spinon lines with disor- 
der, but neglecting the vertex corrections. In the limit 
of low frequencies, n^ 1 + HJ 1 is fully dominated by the 
second term; hence the holons don't affect the residual 
conductivity. We get, on the heavy Fermi liquid side and 
in the limit of low frequencies, 



n 



with n r 



(1 - iu)T c ) 



(44) 



(45) 



In (|44f the effect of impurities is implicitly taken into 
account in the scattering lifetime of the conduction elec- 
trons t" 1 ~ Tq 1 + T 2 . In the limit of zero frequencies we 
get 



m f 



(46) 



while, on the localized side of the transition, only the 
conduction electron conduct. The residual conductivity 
thus jumps at the transition; on the heavy Fermi liquid 
side the f-spinon start abruptly to conduct. This result 
is in agreement with the study o&.. Note however, that 
the jump in the conductivity obtained in this model is 
very unlikely to be detectable , since m c <C m 

We turn next to the temperature dependence of the 
conductivity at the QCP. We first focus on the low tem- 
perature regime, already studied ir*21, saving the study of 
the intermediate temperature regime for the next section. 
As in the previous section, there are two cases of interest. 
In case a) there is a gap in the continuum of the particle- 
hole excitations. Because of the gap, the electron and 
spinon lifetime are not affected by the scattering with 
the holons, and we get the standard Landau Fermi liquid 
law tJ ~ t^ 1 ~ T^ 1 +T 2 . The damping of the holons 
doesn't come from the particle-hole continuum but only 
from the gauge fields This damping of the holons 
itself produces a finite bosonic lifetime. The polarization 
writes (see Appendix [Ej) 



n 6 = — (-iun) 



(47) 



with 



n 
n 



—Log T 



in D = 2; 
in D = 3; 



(48) 



One sees that, in that case as well, the sum H f +H 



is dominated by U h , for D 



2 and by Uj 1 



f 

for D 



Nevertheless, in both cases the conductivity is dominated 
by the conduction electrons as can been seen from (|40[) . 
The resistivity thus varies like T 2 at the QCP. 



p ~ po + T m D = 3 



(49) 



This result contradicts the previous study 24 . In our view 
the contribution of the conduction electrons to the con- 
ductivity was overlooked in 24 . 

The second case of interest is when the two Fermi sur- 
faces intersect, called case b). Here the particle hole 
continuum has no gap, hence the bosons are not perfect 
conductors anymore, but are damped by the particle-hole 
continuum. We write (see Appendix [E] ) 



n b = — (—iun) 



(50) 



where r b carries the temperature dependence of the po- 
larization. We find for z — 2 



U Log T in D = 2; 

n - VT in D = 3; 



(51) 



Now in the Eliashberg theory, damping from the gauge 



fields leads to, 1 



Tq 1 +S/(T) with 

T 2 / 3 for D = 2; 

-TLog T for D = 3 
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One sees that, in that case as well, the sum 14 ^ + Tl b 

is dominated by n^ 1 , for D — 2 and by 14 J 1 for D = 3. 
Nevertheless, in both cases the conductivity is dominated 
by the conduction electrons. Since the regime is char- 
acterized by z = 2, the conduction electrons have the 
inverse scattering lifetime E c ~ J 13 / 2 in D = 3. Since 
the backscattering processes are naturally present in the 
model (see next section) , this regime is equivalent to the 
AF SDW scenario, with 

p ~ po + T 3/2 in D = 3 . (52) 

A last remark for this section is that we recover the 
result of 25 that, although the gauge fields are gapped in 
the Higgs phase, the external electromagnetic field is not, 
which prohibits -thank goodness ! -Meissner effect in the 
heavy Fermi liquid phase. The easiest way to see it is to 
make the change of variable a'^ — a M + in (|36l) . The 
fictitious gauge fields are now a'^ and the external gauge 
fields are decoupled from the holons : 

S mt = 6 f v b • a' b + c) v c • A c + f v f • (a' - A) / . (53) 

One now follows Appendix [F] and write the Ward Iden- 
tity related to the external gauge field. Since the source 
terms can be put to zero ( only the holons get a non van- 
ishing source term) the mass of the electromagnetic field 
vanishes like 

- i<f n A(lA „ = . (54) 



not broken at the phase transition, here translational in- 
variance is naturally broken since the f-electrons are on 
the brink of localization. Hence, there is no need for ex- 
ternal translation invariance breakers, like impurities, to 
break translational invariance. Umklapp processes are 
naturally present which decay the current. In this sense 

Ttr ^ T QP , (57) 

where r tr is the transport time while tqp is the quasi- 
particle lifetime. For z = 3, D = 3, it has the standard 
form 

t qp ~ T ■ (58) 

Our claim is that the electric transport in this phase is 
correctly described by the Drude polarization 




, where here we have used bare lines. More precisely, us- 
ing the Kubo formalism, one can first show that the con- 
duction electron polarization is unaffected by coupling to 
the bosons and the f-fermions. The coupling to f and b is 
protected by gauge invariance. Pictorially we mean that 



X. QUASI-LINEAR RESISTIVITY 



We examine here the transport in the intermediate 
regime. This regime differs from the low temperature 
one by its dynamical exponent which is now z = 3. In 
this regime, the IL composition rules are still valid, with 

U ~ -Log T for D = 2; 
n ~ TVS for D = 3; 

hence the holons still don't participate to the transport. 
As before, the conductivity is dominated by the conduc- 
tion electrons. In this section we examine in more details 
the Drude form assumed in (|44|) . The arguments have 
been given in a previous work, and we recall them here for 
clarity 2 ^. The main idea is that the Drude form is valid, 
with the scattering lifetime of the conduction electrons 
given by 

r" 1 - — TLog T (56) 

in I? = 3. ([55]) is the typical scattering lifetime for a 
z = 3 QCP, like for example a ferromagnet. The im- 
portant question is -"how does the current decay in this 
special q = phase transition ?" The unique feature of 
this QCP is that the current naturally decays through the 
lattice of f- electrons. Contrarily to a usual ferromagnetic 
QCP or nematic QCP where translational invariance is 




+ 



This comes from the fact that [ww] _1 = 14/ + lib. 
The next observation is that the vertex corrections are 
negligible in this regime; namely the next leading order to 
the Drude formulae, is of order (a 1 ) 2 and 1/N, as shown 
in the diagram below. 



15 




[{a'f/N] 



Hence, in this regime the transport is simple and elec- 
trical conductivity can be expressed through the simple 
Drude formulae. The evaluation of the inverse scattering 
time in D = 3 is given in AppendbiHl We get 



Ttr 
Ttr 



TLog(T/E*) mD 
T 2 / 3 in D = 2 . 



(59) 



Note that the logarithm in D = 3 has a thermal origin; it 
differs from the logarithm which appears in the real part 
of the self-energy, in D = 3. 



XI. SUMMARY OF THERMODYNAMIC AND 
TRANSPORT 

To summarize transport and thermodynamics close to 
the QCP, we distinguish the regime I for T < E* and the 
regime II for T > E* . The exponents for transport and 
thermodynamics in regime I depend on the form of the 
spinon Fermi surface. If there is a gap between the spinon 
and electron Fermi surfaces, as is show in Fig[3] for the 
case a), then the anomalous exponents for the effective 
mass - appearing in C v , are due to the massless gauge 
fields with z — 3. The electrical resistivity is dominated 
by the conduction electrons, and since the scattering with 
the spinons is gapped, it follows the usual T 2 law char- 
acteristic of the Landau Fermi liquid. The susceptibility 
doesn't couple directly to the critical modes, hence here 
again the Fermi liquid law is recovered. The exponents 
are summarized in table I. 



T < E* Regime I 


Cv 


Ap(T) 


X(T) 


D = 3 


-TLog(T) 




Xo 


D — 2 


j.2/3 


rj->2 


Xo 



TABLE I: Transport and thermodynamic exponents in the 
low temperature regime when the particle-hole continuum is 
gapped. 



dynamic spin susceptibility, fact which poorly fits the ex- 
perimental observations. 

The second case is when T < E* , but the spinon and 
electron Fermi surfaces intersect. We call this regime I'. 
In that case the particle-hole continuum goes down to 
T = with hot regions at the intersection of the two 
Fermi surfaces. The situation is analogous to the AF 
SDW QCP, except that the critical modes are solely in 
the charge channel. The results are summarized in table 
II. 



T < E* Regime I' 


Cv 


Ap(T) 


X(T) 


D = 3 


-TLog(T) 


rp3/2 


Xo 


D = 2 


rp2fi 


— TLogT 


Xo 



TABLE II: Transport and thermodynamic exponents in the 
low temperature regime, where the spinon and electron Fermi 
surfaces intersect. 



In our view, the most interesting regime is for T > E* . 
In that case, the exponents don't depend on the shape of 
the spinon Fermi surface, but the very existence of this 
regime requires the presence of a spinon Fermi surface 
at the QCP. Here the transport is simpler than closer to 
the QCP. The resistivity shows a quasi linear behavior in 
D = 3 and a sub- linear exponent in D = 2. The results 
are summarized in table III. 



T > E* Regime II 


Cv 


Ap(T) 


X(T) 


D = 3 


TLog(T/£*) 


TLog{T/E*) 


T 4/3 


D = 2 


rp2/3 


T 2/3 


-TLog(T) 



TABLE III: Transport and thermodynamic exponents in the 
Maginal Fermi liquid regime around the KB-QCP. The expo- 
nents are in agreement with those of Ref.— . 



Note that the temperature dependence of the spin sus- 
ceptibility, although departing from the standard Landau 
Fermi liquid law, still doesn't show anomalous exponents. 



XII. CONCLUSIONS 



We see that the situation is peculiar in the sense that 
the Landau Fermi liquid paradigm is broken for the ther- 
modynamics only. The resistivity behaves as T 2 even 
though the residual resistivity jumps at the Fermi sur- 
face. There is no trace of anomalous exponents for the 



In this paper we have given the simplest consistent 
treatment of a selective Mott transition in the Ander- 
son lattice, using a U(l) slave boson technique associ- 
ated with an Eliashberg treatment of the vertices. We 
find that the QCP exists and has a multi-scale charac- 
ter. At the QCP below E* , the exponents for thermody- 
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namics and transport depend on the shape of the spinon 
and electron Fermi surfaces. When the particle-hole con- 
tinuum is gapped, the anomalous transport scattering 
is gapped as well, and the resistivity follows the landau 
Fermi liquid law. The effective mass is dominated by 
the fluctuations of the transverse gauge fields showing a 
Reizer singularity^. If the spinon and electron Fermi 
surfaces intersect, the particle-hole continuum is gapless 
and the transport and thermodynamics are anomalous 
down to the lowest temperature, with z = 2 for T < E* . 
Above the energy scale E* , the resistivity doesn't de- 
pend on the shape of the spinon Fermi surface, and we 
get a universal quasi- linear resistivity in D = 3. In our 
view, the important question raised by our study is the 
question of the presence or not of a selective Mott tran- 
sition in the Anderson lattice. Said in general words, one 
can reasonably ask whether we believe the anomalous 
properties observed in heavy fermions are due to a Mott 
localization of the f-electrons. If we believe it is the right 
answer, then the U(l) slave boson gauge theory is the 
most straightforward approach to such a phenomenon. 
It would be very interesting to have studies from other 
techniques, like DMFT, to first confirm the presence of 
the transition and, if it there, to give more details about 
the elementary excitations. 

It is reasonable to ask whether the U(l) slave boson 
theory, although being the simplest description of the 
Mott transition, is the appropriated tool. This question 
is of relevance as well for the cuprate superconductors, 
where gauge theories, with sometimes bigger algebra like 
SU(2), have been extensively used to describe the ap- 
proach to the Mott stat o i 45 . It is clear from the above 
study, that this approach suffers from what we would 
call " spinology" . Fermionic spinons with a finite Fermi 
surface are needed at the QCP for the QCP to exist at 
all. One can reasonably question this feature, and won- 
der whether under more powerful techniques, this fea- 
ture would survive. Nevertheless, before discarding the 
U(l) slave boson gauge theory for the Anderson lattice, 
one must consider that it gives a very interesting regime 
with a quasi-linear scattering and transport lifetime in 
3D. This unique feature is not easily obtained within 
any theory and the good point is that it doesn't depend 
on the shape of the spinon Fermi surface, but only on 
its presence at the QCP-hence it is a direct consequence 
of the "fractionalization" of the heavy electron at the 
QCP. We can use this regime for making experimentally 
testable predictions. A first application to He 3 bi-layers 
has been performed 31 . We can as well make predictions 
for the thermal transport in this regime, and call for ex- 
perimental confirmation or invalidation. 

I am indebted to M. Norman and I. Paul for count- 
less discussions on the similar Kondo-Heisenberg model. 
Many ideas of this paper came from interactions with 
them, related to this previous study. The description of 
the quasi-linear resistivity benefitted from a useful dis- 
cussion with D.L. Maslov and the study of the Ward Iden- 



tities and the Ioffe-Larkin composition rule, from useful 
interactions with K-S. Kim and O. Parcollet, and a very 
useful, heated discussion with P. Coleman. A special 
thank to C. Bena and A.V . Chubukov for a useful read- 
ing of the manuscript. This work is supported by the 
French National Grant ANR36ECCEZZZ. 
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APPENDIX A: EVALUATION OF THE 
INTEGRALS FOR THE MEAN-FIELD 

At T = 0, the calculation of the integrals used in the 
mean- field treatment is analytical for linearized bands: 
We call 



A 
B 
C 



T ^ G ff( k > iuJ n) 

fe,(7,U> n 

T ]T G fc {k,iu> n )/{bV + a ) 

k(T : LJ n 

T ^2 e k G ff (k,iuj n ) . 



We diagonalize the 2x2 matrisi^ which accounts for the 
hybridization of the f- and c- bands: 



Xo 



4 + eu 



e° k + e k + VA 



A = {e a k -e k f+A{a Q + bV) 2 . 

The integrals are all performed in the same way, first by 
summing over the Matsubara frequencies , and second 
by doing the momentum integration. The momentum 
integration is done by linearization of the band. 



A = NT^, 



(iuj n - e k ) 



(iu> n - Xi)(iu> n - X 2 ) 



Np Q I de 

-D 



~n F (z) 



(z-e) 



dz 



2in (z - Xi)(z - X 2 ) 
where the contour is on the whole complex plane 
" fn F (X 1 )(X 1 -e k ) n F (X 2 )(X 2 - e k ) 



Npo 
2 j-u 



-d V (Xi - x 2 ) 



(X 1 - X 2 ) 



de 



de- 



-y + ^y 2 + ±{bv + <jy) 2 

0/ 2 + MbV + a ) 2 
-y-^V 2 +A{bV + a f 
Vy 2 + MbV + *o) 2 



with e m and e p the Fermi levels for the upper and lower 
bands respectively. 



e m = (-£/ + a'p - yj ( £f + a'p) 2 + Aa> '(bV + a a ) 2 )/(2a') 

e p = (-£/ + a'p+ ^j{e f + a' p) 2 + Aa>{bV + a a ) 2 ) / (2a') 

with the conditions — D < e m < 0; < e p < D , 
and a' = ab 2 + 4>o/D . 

One obtains 



A 



(~2y^ D + y m - ^y 2 m + A(bV + a«f 
1 — a') \ 



2(1 - a 1 ) 



+ y P + ^y 2 + <bv + a ) 2 ) , 



Vm = (1 - a')e m - e f - p; 
where y p = (1 - a ')e p - e f - p; 

y-D = -(1 - a')D -ef- n . 

We proceed in the same way for B and C to find: 
Np 



B = 



(l-a>) 



7S L °9 



y P + ^y 2 P + ^bV + a ) 2 



D 



2(bV + a ) 2 Log 



+ ( e f + v)ym + yl 



y m + Vyl + HbV + a^ 
y P + yjy% + *0>v + °o) 2 



- (y m /2 + e f + pWyl + i(bV + a Q ) 
+ (e/ + M)y P + Vp 

- (y p /2 + e f + p)^y 2 + 4(bV + <j ) 2 



APPENDIX B: EVALUATION OF THE 
INTEGRALS FOR THE VERTICES 

In this section, we focus on the evaluation of the ver- 
tices for the regime (ii) where we have logarithmic fre- 
quency dependence of the polarization. We call g — V 2 po 
with V the coupling constant between the soft modes and 
the f- and c-electrons and po the density of states of the 
c-electrons. Note that although it is not possible to form 
the vertex correction at the fist loop, we still evaluate the 
fictitious one-loop vertex below, knowing that the two- 
loops vertices behave as products of the one-loop ones. 
We start with the static vertex: 



r(o,o) 




-NLog(a') J (27r) d+1 q 2 - aLog\u n \ 
1 



(iuj n + a'v F qcos9)(iu> n — vpqcosQ) ' 

g f q 2 dqdu! 1 



NLog(a')(l - a') J (2tt) 3 -iu> n (q 2 + aLog\w n \) 

x = l 



— ) Log 

VFq 



v F q — x 



a' VFq — x 
g(-Loga) 



x=-l 



N(l - a')l2(2Tr) 2 E F 



18 



with a = 1/ Log(a'). Hence the static vertex is small in 
1/N. 

We next evaluate the dynamical vertex, with linearized 
Fermi surfaces 



r(q,n) 




r(«, fin) - _ NLog ^ J ( 27r )d+i q 2 - aLog\cj n 
1 

(iu> n + a'vFqcos9)(iuj n + iCl n — vpqcosQ — vfQx) ' 

-NLog(a')(2ir) 3 iCl n - vfQx 7_n„ Log\u> n \ 
9 



-NLog{a'){2irf 



ifi n - v F Q x 



(Q, n LogLog\Q, n \ - Li(—Q n )) 



We see from ()B1|) that for Q x — the vertex has a 
LogLog- singularity. It is not small. 

The same evaluation with the curvature of the Fermi 
surface gives: 



.9 



d d qdu! 



1 



-NLog{a') J {2n) d+1 q 2 - aLog\iv n \ 
1 



(iu> n — a! vpqcosQ — a! q\ / (2m)) 
1 

(iuj n + ifl n — vpqcosQ — vpQ x — q 2 ^/(2m)) 
(—iduj) x 



9 

-NLog(a')2(2ir)* J_ n J 
1 



-Logo 1 



{QxVf + (1 - a')Log\uj n \/ 'Loga') ' \Log\uj r , 

9 itt n Log(~Log\n n \) 

-NLog{a')2{2Tr) 3 (1 - a' Log\Q n \/ Loga') ' 



We see now that the curvature regularizes the vertex both 
in large N and in the infra-red frequency sector. 



APPENDIX C: INSTABILITIES BEYOND THE 
ELIASHBERG THEORY 

In this section, we evaluate diagrams beyond Eliash- 
berg theory, but that are potentially dangerous for the 
static sector of any q = QCP. As mentioned in the 
main text, such singularities were discovered by BKV 18 
and are typical of the type of problems coming from the 
presence of a finite Fermi surface in the theory. Precisely 
we want to evaluate 



n£(T) 



n, 



(i) 




We first note that the two diagrams are proportional: 
Ili 1 ^ = a' and that there is no corresponding vertex 
insertion at the first order. One can check that, in the 
low energy regime T < E* , the two diagrams are not 
singular, since the average gap between the spinon and 
electron Fermi surfaces protects it. Here we want to check 
the stability in the intermediate regime for T > E* . We 
have z = 3 in the boson propagator. To understand the 
source of the problem it is instructive to compare the 
following four-field diagrams 



/ T 9F i f 



and 



9i 

X 



h 



where g± is a mode- mode coupling constant, coming for 
example from the term —J riinj /4 in Eqn ([3]). The 
54 mode-mode coupling is standard in a 4 -theory and 
provides correction to scaling extensively studied in, for 
example 62 chapter 42. If 34 > the </> 4 theory is stable 
and, close to a QCP, the corrections to scaling follow the 
law 



TOb(T) 



b b 



rp(d+Z-2)/z 



One sees that, in a fermionic theory, one can form correc- 
tions to scaling from the fermion vertex gp which leads to 
our two diagrams ni 1 -* and . The difference between 
gp and 34 is that the fermion loop is dangerous, and can 
change the sign of the vertex. According to the value of 
z, it can as well lead to a more relevant term than the 
standard 4> 4 corrections to scaling. 

We turn to the computation of the diagrams. With 
g = 8k F V 2 /p , c = p N/(a'v F ). 

U^(T)=NgpoT I 4zL d ^D b (q,{l m ) 

x G 2 c (k,u n )Gf{k,w n )Gf(k + q,w n + fi m ) ; 



q + aq z 



(iuj n + i£ c (w„) - e k + fi) 2 (iuj n + ffi/(w„) - a'e k - e f ) 
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(iu n + iCl m + iE/(w„) + iS/(f2 m ) - a'e k - a'v F q x - £/) 

We have Sgn[E/ )C (a; n )] = Sgn[aj„].We perform first the 
integration over , noticing that the integral is non zero 
if the poles are in the same half plane. To fix the ideas 
we take ui n < while uj n + Cl m > (the result is identi- 
cal with the other choice). We then close the contour in 
the upper half plane to catch the pole coming from the 
last factor. We obtain, after neglecting the terms propor- 
tional to a' and neglecting iu> n compared to tEf tC (ui n ) in 
the Green's functions: 



Jq 2 x + q 2 y ^\q y \+ql/(2\q y \). Wc find 



H^(T)^Na'gp T 

2vk 

(ffi/(f2„) - a'v F q x ) 
1 



d d q 



1 



(iE/(fi„) + ffl/(w n ) - a'v F q x - e/) 2 



Na'gpoT 



d d q iO m 
{2n) d c\n m \/q + aq 2 



(CI) 



(iE/(f2 m ) - a'v F q x ) 3 
The calculation in D = 2 and D = 3 differs at this point 



n«(T)~-r 4 / 3 



(C2) 



m/0 



I?* I 



x / rffe (|9„|+d/(2|«vl)) 



c|0„ 
1 



(iE/(f2 m ) - a'v F q x ) 3 



the integration over g y now leads to a logarithmic singu- 
larity in q x 



1 sing 



= I dq x 



ql Log(A/| ga |) 
(iE/(Q m ) - a'v F q x ) 3 



Keeping only the dependence in S7„ in the integrand, we 
get 



— fsing is performed by continuation in the upper half 
a ^plane, if Q m < and in the lower half plane if £l m > so 

that to avoid the pole in the Green's function. Changing 

variables in q x = iz we get 

Ising = - / ldzSgn{n m )— —————— -3- : 

Jo (-i) 6 (|E/(S2 m )| + a'v F z) 

-i7rSgn(» m ) / A 
{a'v F f H|S/(fi ro )| 

Putting things together we get 



D 



Since there is only one pole in the last factor, the inte- 
gration over q x lead to a typical q x of the order of Ey (it 
is the reason why full fermion lines can be safely replaced 
by bare ones in the Eliashberg theory). One possibility 
is that the integration over q x doesn't vanish due to the 
branch cut in the boson propagator q x ~ q± ~ 117ml 1 / 3 . 
It is what happens in D = 3. We find 



Note that the minus sign is of importance. We argue 
though, that since we are in the intermediate regime 
where T 3> E* , the correction in T 4 / 3 is irrelevant, just 
giving an extra UV cut-off for the intermediate regime. 
The stability of this regime is thus ensured in D = 3. 



D 



In D = 2 the situation is more complex. The branch 
cut at q x ~ q ± ~ \Q m \^ 3 would lead to U { a 1} {T) ~ — T, 
but there is a stronger singularity first discovered in 
Refi 1 ^. Indeed, we suppose q y q x and expand q — 



n d)(T) = n^^tY 

V* (2tt) 2 ^ c \n 

y 1 m^O 1 



-z7rSgn(Sl m ) 
(a'v F ) 3 



Log 



A 



|S/(fi m )| 



since TJ^^a/t = we nnc ^ that the T-dependence 
of the above comes from the m = term only. Finally 



n«(T) = 



q'gpo / A 

"(2tt) 2 c(a'v F ) 3 



(C3) 



This result puts the intermediate regime in a fragile sit- 
uation. This term is of negative sign and dominant com- 
pare to E*\ it has the potential to destabilize the regime. 
Note that in D — 2, corrections to scaling coming from 
the 54 interactions have the same temperature depen- 
dence (—TlogT), but with a positive sign. Corrections 
to scaling hence compete with Il£ 1 ' > and tend to stabi- 
lize the fixed point. We have reached here the limit of 
the Eliashberg theory. Whether the intermediate regime 
is stable or not in D = 2 depends on whether we have 
strong enough corrections to scaling of positive sign. This 
requires, for example, strong enough ferromagnetic short 
range fluctuations (Jo < in Eqn.Q). The stability 
of the intermediate regime is a matter of pre-factors be- 
tween the two terms. 



APPENDIX D: IL COMPOSITION RULE FROM 
THE CONSTRAINT 

We can recover simply the IL composition rules by ap- 
plying the constraint nf + rib = 1. The external current is 
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associated to the conduction electrons and to the bosons 
b. 



J — J c + Jfc . 



(Dl) 



but introducing now the fictitious field e and the external 
field E we get 



J c = a c E 

J fc = <7b (E + e) 

J f = o } e . 

We apply J/ + J b = to |D2j) and get 



E 



Of + Ob 



which leads to 



O = Or 



Of + <Jb 



(D2) 



(D3) 



(D4) 



APPENDIX E: EVALUATION OF THE 
POLARIZATION AND SELF-ENERGY OF THE 
BOSONS AT THE QCP 

1. Self-energy 



£&(fi n ) 



This self-energy captures the effect of the gauge fields 
on the boson propagator. This effect is sub-dominant 
in all regimes but in regimes I where, because of the 
gap in the particle-hole continuum, the only source of 
damping for the bosons are the gauge fields. We evalu- 
ate Ef, in regime I where T < E* and the particle-hole 



continuum is gapped. Renaming c — > Nir/(2r 
a — > l/(2m/fc|i), we have 

d d q 



l f v f. 



and 



(27r) d 2dm 2 —iu n + aq 2 



c\uj„ + n n \/q + aq 2 



(El) 



Performing the analytical continuation and integrating 
over the two branch-cuts uj n = and u) n + Q n = 0, we 
get for fi„ > 

U) \ d d q q 2 



/^Cothf^) 

J tilt \2TJ (2ir) d 



2dm\ 



1 



1 



(— uj + aq 2 — id) (c(—iuj + fl n )/q + aq 2 ) 
1 1 



(— u> + aq 2 + iS) (c{iui — Q n ) / q + aq 2 ) ) ' 



*L coth (±L\ J. 



an — V2t; j^2d^ b ^ )5{ ~ UJ+aq2) 



Coth 
2aq 2 



a 2 q 4 — [icujq — cfl n /q) 2 
-1 r d d q 
~2~ J (2ir) d 2dm 2 



q2 Coth faq2 



2aq 2 



% \2T J (icaq- cfl n /q) 2 

Since this integral is dominated by large q, we get 

E h (n ) ~ -\n |( d + 2 )/ 2 [ A — - — 

b{ n) ~ 2 1 n| Jv2T(^) d Mm 2 (ix 2 -l) 2 

= (A + iB) |n„jC d + 2 5/2 1 (E2) 



where 




r A d d x a x e 
JO (2n) d idm% (ix'^-iy 

f A d d x a x 6 
JO {2ir) d idmi (ix 2 -l) 2 



in the limit where T — > 0. We use in the next section the 
short-hand notation: 

£ 6 (O„)=/ a |O n r with " q = _ (d , + L 2) ^ 2; ■ (E3) 

Jq — A + in 



2. Polarization 

In the following, we evaluate the polarization bubble 
of the bosons in the three possible regimes at the QCP. 




II&(g,ifi n ) 



The generic form of the bosonic polarization is 

n b (Q,zao = T^ j -0 



2 2 



(27r) d 

-D fc (q, iuj n )D b (q + Q, iu n + iQ n ) , (E4) 



with vt, the vertex defined in (|36[) . 



3. T = 0: form used in the gauge propagator 

We evaluate here the bosonic polarization contributing 
the the gauge fields propagator. Here Vb — l/m&. At 
T = only undamped holons contribute to the damping 
of the polarization 



1 



1 



d d q q 2 
(2ir) d 2dm\ -ioj + aq 2 



-(iu> + iQ n ) + aq 2 + qQCosO /mb ' 
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with a = l/(2m&) 

2iir f Vld-iq d dqdCos6 q 2 — O n + iqQCosO/mf, 



-i J (27r) d 2dm 2 b ft 2 + (qQCos9/m b ) 2 

where Qd-i is the solid angle of dimension d-1 
il^q^dq q 2 7r|Q n | 



(27T) 



<i-l 



2dm;, Q 



4. form used in transport 

gime i/iere is a gap between the two Fermi surfaces 



The form of the boson propagator is given by (J34J) with 
Ufa given by the frequency dependence of (|19|) and 
given by (|E3[) ■ We see that the particle hole contribution 
IT/ C to the polarization is undamped in this regime; the 
only source of damping £& comes from the gauge fields. 
For the Kubo formulae, we evaluate the polarization at 
q — ► 0, we need only to retain the damping part of the 
polarization Y, b . The boson polarization then writes ( 
with c = 1 and a — l/(2m&), v b — l/m b ) 

n fc (0,ztt) = r£j Wf 2dm 2 f^ n \»+aq 2 
1 



f^\u n + n n \ a + aq 21 

considering the two branch cuts at u> n = 
and uj n + fl n — we get, 



f duj _ , / cj \ d d g 

/ Coth — - — f 

7 4i7r V2T/ 2tt ' 



x 



27V (2?r) d 2dm 2 

V/ tt (-^)« + aq 2 ) if$(-iu> + 0„) a + ag 2 ) 
1 1 



/ Q M«+ag 2 ) 


(/o Q (-^ 


+ o„)« 


+ ag 2 ) 


1 




i 




ftf(-iuj) a +aq 2 


) (/o a (™ 




+ ag 2 ) 


1 




1 


aq 2 )) 


/ Q M« + ag 2 ) 


(fo(^ + 


ao Q + 





with the tadpole diagram). We get 

to \ d d q q 2 



U b (Q,iQ n ) ~ / ^- Cothf 



4i?r V277 (27r) d 2dm 6 

2iSin[air/2] ftfuj a 



2 



((/ Q w Q Cos [an/2] + ag 2 ) 2 + (Sin[an/2]f*uj a ) 2 ) 3 
x [-2w Q - 1 Cos[(a - 1)tt/2] ((f£uj a Cos[air /2] + ag 2 ) 2 

- (Sm[mr/2]/o a ^) 2 ) 

- 2iw Q " 1 S'm[(a - 1)tt/2] (-iS'm[a7r/2]/ Q w Q ) 
x (/ Q w a Cos[a7r/2] + ag 2 )] . 

From the above formulae and using a = (d + 2)/2, we 
extract r b in all dimensions to get 



Tb ~ —logT in £) = 2 
t 6 - T 5 / 4 in L> = 3 



(E6) 



b. Regime I : the two Fermi surfaces intersect each other 

We start now from (|24|) (where we have re-named 
Pqc/ {a'vFqo) —> c and poa — > a, v& = l/m b ). We take 
fi„ > 0. 



n 6 (o,in„) = T^/A^ 



d d g g 2 



1 



(27r) d 2dm 2 c\to n \ + aq 2 



1 



c|w„ + n n | + aq 2 ' 

considering the two branch cuts at u> n = 
and uj n + Q n = we have , 



= /" C,nt h (JiL} -^i, 

1 



Coth 

1 



4i7r V27V (27r) d 2dm 2 



(— ciw + aq 2 ) (— c(iw — Q n ) + aq 2 ) 
1 1 



(ciuJ + aq 2 ) (—c(iu) 


-n n ) 


+ ag 2 ) 


1 


i 




(— ciuj + aq 2 ) (c(iu> 


+ Q n ) 


+ ag 2 ) 


1 


i 


- «9 2 ) ) 


(ciw + ag 2 ) (c(mj + 


On) 4 





<^ p d d q q 2 

4nr \2TJ (2n) d 2dm 2 



1 



1 



f°(-iLu)<* + aq 2 /«M« + ao 2 

1 1 



f«(-iLU + n n ) a + aq 2 f£(iuj + n n ) a +aq< 



We expand the second factor in fi n and take the part 
proportional to |f2 n | ( since the constant part cancels 



coth ^i- 



4wr ~ "\2TJ (2ir) d 2dm 2 c 2 u 2 + a 2 q* 
2cO„ + ag 2 



c 2 w 2 + (cO„ + ag 2 ) 2 ' 

Since the bosonic mass cancel out with the tadpole dia- 
gram (see Appendix [Fl or [G|l . we have to extract the part 
proportional to f2„, to get 

duj / u> \ d d q q 2 c 2 LoQ n 

~ \2TJ (2^p 2dm| (c 2 u 2 + a 2 g 4 ) 2 ' 
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From the above formulae, we extract Tf, in all dimensions 
to get 



Tf, ~ —logT in D = 2 
n ~ T 1 / 2 in D = 3 



c. Regime II both cases 



(E7) 



In this regime we start with Eqn. p2|) for the bo- 
son propagator, showing z = 3 (with the renaming 
Po/(o:'vf)c — ► c and p^a —> a) . We take as well f2„ > 0. 
Here Vb is a bot more comple since we have to expand 
|w„|/(|<7 + a|) in the first order in the vector field a to find 
Vb- We get 



Vb = l/m b - cuj n /q 3 
Since in the integral below 



(E8) 



n 6 (o ) ffl n )=T^ J 



d d q q 2 v t 



2„,2 

b 



(2ir) d 2d c\u n \/q + aq 2 



1 



c\uj n + Q n \/q + aq 2 ' 

considering the two branch cuts at u> n = 
and oj n + fi„ = , 



2, ,2 



4ct V2T/ (2^) d 2d 

1 f 



(—cicu/q + aq 2 ) (—c(iuj ~ Q n )/q + aq 2 ) 
1 1 



+ 



(ciuj/q + aq 2 ) (— c(icj - 


-fin)/tf- 


- aq 2 ) 


1 


1 




(—ciuj/q + aq 2 ) (c(iu - 




-aa 2 ) 


1 


1 





(ciuj/q + aq 2 ) (c(iuj + Sl n ) / q + aq 2 ) _ 



APPENDIX F: WARD IDENTITIES (WI) 

In this section wc derive the WI associated to the gauge 
invariance of our theory. When it is not mentioned, the 
field a has been put to zero at the QCP. The first point of 
interest in to show that the mass of the gauge fields is zero 
in the Coulomb phase. The second point is to show that 
it is non zero in the Higgs phase. In the case of a purely 
bosonic gauge theory, this second point is simple and can 
be found, for example, in the Peshkin-Schroder 52 . In the 
case of gauge theory with non relativistic fcrmions, the 
point is more subtle since the gauge fields a M , a v are not 
only coupled to the Higgs boson. Gorkov 61 was the first 
to show that the mass is generated in the Higgs phase, 
in the case of a supercondcutor. Here we follow Zinn 
Justin on page 432 62 , with a field theoretic derivation of 
the same result. The possibility of massless Higgs phase, 
although non generic, will appear at the end. 

We start from Eqn. (f2"T)) with the gauge fields both cou- 
pled to the spinons f and holons b. For simplicity we have 
put a = everywhere since this parameter is irrelevant 
at the QCP. 



So = 



- J d d xdr 
fl(x) (dr + 



+ X + E )U(x) 



(V - iea/c) 2 
2m f 

+ 6t(l) ( aT+ (v^ + ,) 6(l) 

+ / d d xdr{b{x)fi{x)c{x) + h.c.) + H c . (Fl) 



For each field, we introduce a source-term, like 
a^i ► J , 

/-»?; 
fi ^ 

b -> J b ; 
tf^Jb 

so that 



(F2) 



c th f—\ ddq q2v t* 2ciuj /q 

Aiir \2TJ {2n) d 2d {c/q) 2 LU 2 + a 2 q 4 

2cQ n /q + aq 2 



{c/q) 2 uu 2 + (cn n /q + aq 2 ) 2 ' 

and taking the part proportional to Q n , we get 

—C th (—\ q2v * (c/q) 2 ujfl n 

tf ° \2T) (2Tr) d 2d {{c/q) 2 io 2 +a 2 q 4 ) 2 ' 

From the above formulae, we extract Tf, in all dimensions 
to get 



T~b ~ ~logT in D = 2 
Tb-T 1 / 3 in L> = 3 



(E9) 



S — So + S s 



, with 



(F3) 



d d xdr (a^Jf, + fir) + fjf + J b b + J b ) 



The part of the action So ()F1|) is invariant under the 
gauge transformation (|2l)|) . Only the source terms S source 
are affected by the gauge transformation. Using the lin- 
earized form of the U(l) algebra e 10 = 1 — iff we get for 
the variation of the source terms: c 



ss. 



source 



d d xdr 



^ + ie(fjf-fiv) 

e up 



i9(bJ b -tfJb)) 
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Now from the change of variables 

&->6(l-i0), 
a,, — > a„ 



(F4) 



;>j:,. 



we check that the whole action S = Sq + S source is 
invariant under the U(l) gauge transformation. Hence 
<5 S source — 0. Integrating by parts the first term in (|F4|) . 
one gets one version of the WI: 



d^ + ieifjf-p^+ieiJbb-tfjb) Q =0 



(F5) 



(|F5j) is applied to any generating functional for correla- 
tion functions. We can first use the generating functional 
of the source currents W{J, rf). 

n /_<9IF dW \ fdW - <9IF\ „ 

(F6) 

But we can also use the generating functional for the 
vertices T, where T is a p-leg vertex, which leads to the 
more useful WI: 



(F7) 

To get a better idea, let's derive the WI for the two- legs 
vertex, which is nothing but the total polarization W^ v . 
We differentiate (|F7[) with respect to a v {y) to get 



-dp. 



d 2 T 



da v (y)da^{x) 
d 2 T 



da v (y)df(x] 

d 2 T 
da v {y)db{x) 



b(x)-tf(x) 



d 2 T 



da v {y)dp{x) 
d 2 T 



da u (y)dtf(x) 



(F8) 



. 



To get the proper vertices, we Fourier transform and then 
put the sources to zero. We note here that it is possible 
to put the sources to zero in the Coulomb phase, but not 
in the Higgs phase where the boson acquires a non zero 
mean-field value. There are two cases of interest : (i) the 
Coulomb phase and (ii) the Higgs phase. We get (i) in 
the Coulomb phase : 



-iq»T 2 (q,-q) = 



(F9) 



which is re- written, with the notations of the body of this 
paper 



iq»Tl flu {q,-q) = 



(F10) 



From (|F10|) we see that, in the Coulomb phase, the WI 
ensures that the mass of the transverse gauge field prop- 
agator is zero to all orders. Note that the identity (|F10[) 
constraints only the mass of the vector fields; since the 
mass is taken at q° = uj = 0, for which the scalar field ay 



is dropping out of (|F10[) . Now (ii), in the Higgs phase, 
the WI writes 

- iq^l^q, -q) + ie(b) (T 2 bl/ (q, -q) - T 2 fl/ (q, -q)) = 0. 

(Fll) 

We see that the gauge field propagator gets massive as 
soon as r^(q, —q) — T 2 fi/ (q, —q) ^ 0. This phenomenon 
is nothing but the Meissner effect for superconductors. 
Note that it can happen that, for special forms of the 
Fermi surface, T 2 v {q, —q) — T? fii (q, —q) = . We then get 
the equivalent of massless superconductivity, for a U(l) 
gauge theory. 

We can also derive the same kind of WI for the three 
legs vertex. Let's take for example 



_«fy(g)y fk 

Jk+q 

We differentiate (|F7j) with respect to f(y) and then to 
f'(z) and put the source terms to zero, except in the 
Higgs phase: 

-d u „ , (F12) 



da,(x)dp(z)df(y) 
d 2 T 



d 2 T 



df(z)df(x) xy df(y)dp(x) 
ie b(x)- 



dff(z)df(y)ab(x) 



= 



dp(z)df(y)dtf(x) 
(i) In the Coulomb phase 

-iq^ fl ^ fk +ie(G- 1 (k + q)-Gj 1 (k^ = 0. (F13) 
In the Higgs phase (ii) it comes 

-^fi + J k +^(G 1 \k + q)-G-\k)) 
Q3 T Q3 T 



3(6) 



db q dft +q dfk db\dftdf k+q/ 



= 



Quite generically, the p-legs vertex is related to the (p-1)- 
legs vertex through the WI. Note that a relation similar 
to (|F13|) can be established for the ferromagnetic QCP, 
using the translational invariance instead of the gauge in- 
variance; one can follow the same steps using the Noether 
theorem associated with translation invariance, instead 
of the equivalent of it, associated with U(l) local invari- 
ance, that we derived here. 



APPENDIX G: DIRECT CHECK OF THE 
VANISHING OF THE MASSES 

In this section we directly check that the masses of the 
gauge field propagator vanish in the Coulomb phase, at 
the first order in perturbation theory. 
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1. Fermion part 

Let's start with the fermions and check the following 
cancellation: 

/ 



11/(0,0)= w» 

f 



n /(°.°) = ^E/(fc 



a u a,. 



d d k k 2 1 



) d d (iu> n - 



2m/ ^ J {2ir) d (iuj n - e° k ) ' 

with e° = fc 2 /(2m/) + e/. Relating the hrst term to the 
second through the identity 



(-2m 2 /) 



d 1 



9m/ (iu n - eg) (iuj n - e° k ) 



we get n/(0,0) 



Noticing that 



1 d 1 



ddrrif 2m / / 7 (27r) 



d d k 



(Gl) 



d d fc 



d/2 



we finally get at the first order in perturbation theory 
11/(0,0) =0. (G2) 

2. Holon part 

The bosonic part follows the same steps as the 
fermionic one. 



n 6 (o,o) 




a u a,. 



A q 2 



1 



, ^ ... ,_„/ d (tw„ - <? 2 /(2m b )) 2 

T ^ f d d q 1 
+ 2^~ b j^J (2ir) d (ioj n - q 2 /{2m b )) ' 

Relating the first term to the second term through the 
identity 



d 



1 



we get lit (0,0) = - 



1 d 



1 



d dm b 2m b 

2 



Noticing that 



d q 

^-^n B (q 2 /(2m b )) 



A± [ n b {q 2 /(2m b )) 



with Qd the solid angle in d -dimensions, we finally get 
at the first order in perturbation theory 



n 6 (o,o) = o 



(G3) 



APPENDIX H: CONDUCTION ELECTRON'S 
LIFETIME 

In this section we give the evaluation of the tempera- 
ture dependence of the imaginary part of ImE c . We work 
in D = 3, and 



d d q q 
(2ir) d c\ui n — ttT\ + aq 3 iw v 



1 



z k+q 



(HI) 

where the frequency ui n = (2n + \)itT. We close the 
contour around the two branch cuts for uj n = and uj n — 
ttT = to get 

E-CT) = -9 I— Cothf-i) /-^ 
V ; y J Aiir \2TJ J (2ir) d 



1 



1 



—icuj + aq 3 uj + i7rT — ai>Fq. 
1 1 



iccj + aq 3 uj + i7rT — avFq x 







f — Tanh 


k2t) J 


/ Aiir y 





1 



1 



icuj + c7rT + aq 3 uj + i5 — avFqx 
1 1 



icuj + cirT + aq 3 uj — i5 — avpqx 



doj f uj 

-g I — — Coth 



4l7T 

2icuj 



\ I qd d q 
2TJ J {2TT) d 
1 



c 2 ui 2 + a 2 q b uj + iirT — avFq x 
^Tanh^ f**2- 



4/;t V277 J (2?r) d 

2 1 
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Taking the integration over q x leads to 



Sc(T) 



avF 



OtVp 



^Cothm 
4tt V2T/ 



An 



Tanh 



V2T/ 



Q q q( d ^dq icuj 
(2 7 r)( rf - 1 ) c 2 w 2 + aV 
n q q^dq 1 



where fid is the solid angle of dimension d. This integral 
is dominated by the low energy part of the first term (the 
high energy part of the first and second terms cancel out) 
which leads to 



ImS c (T) = 



.9 



avF 



T duT 



where am is a IR cut-off. Taking d 
variables for x = q 3 /oj we get 



( 27r )(d-l) C 2 W 2 +a 2 q 6 ' 

3 and changing 



Im£ c (T) 



fJ 



avFC 
TLog 



am 

(- 

\a-IR 



T r°° Andx 1 
W wi 3(2^)3 l + a 2 x 2 



(H2) 



The question is now to determine the cut-off am. Since 
we work at finite temperature, there are two sources of 
IR cut off which are E* and irib(T) , with rrib(T) is the 
temperature dependence of the holon mass. 

a IR = Max[E* ,m b (T)] . 

TOf,(T) is determined by evaluating the corrections to 
scaling 



b 




where is the coupling constant associated to the tfi 4 - 
holon field theory, if it is there. One can check that the 
first diagram goes like T^ d + z -^l z the second one like 

T (d+2)/2 and thc third Qnc Hke T 5/3 m D = 3 Hence in 

D = 3 

m b (T) ~ T 4 / 3 . 

Note that the lines are full, and that we must use the 
Eliashberg theory for this check. Hence in this model, in 
the intermediate energy regime, TOb(T) <C E* . We get 



L null's r, 7'L,.,,( Jl 



(H3) 
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